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ABSTRACT 

The present paper describes MHD free convection flow past inclined porous stretching sheet under the 
influence of viscous dissipation and Radiation effects. Similarity transformations are used to reduce non linear partial 
differential equations to ordinary differential equations. Approximate solutions have been derived for velocity, 
temperature and concentration, using R-k method of fourth order along with shooting technique. Graphical analysis has 
been done to identify influences of different physical parameters on velocity, temperature and concentration. 
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1. NOMENCLATURE 

OC Angle of inclination 

5q Magnetic field strength 

<7 Electrical conductivity 

Gy Stefan-Boltzmann constant 

p Density of the fluid 

Xj/ Stream function 

U Velocity component along x direction 

V Velocity component along y direction 

T Fluid temperature 

6 Dimensionless fluid temperature 

7! Dimensionless distance normal to the sheet 

7^ Free stream temperature 
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T w Uniform sheet temperature 

K Thermal conductivity 

Cp Specific heat at constant pressure 

Qq Volumetric rate of heat generation/absorption 

q r Heat flux 

C Concentration 

Free stream concentration 
(p Dimensionless concentration 

C w Concentration of the sheet at wall 

K Thermal conductivity 

/fj Mean absorption coefficient 

V Kinematic viscosity 

g q Acceleration due to gravity 

k Porosity 

2. INTRODUCTION 

Investigation of radiation interaction with natural convection got good significance in recent years. This is due to 
the presence of thermal radiation in the surface heat transfer. In technical processes like paper production, crystal growing, 
glass blowing, aerodynamics, we can observe boundary layer flow on continuous surfaces. Because of this, heat transfer 
analysis plays vital role in industries. Ostrach [9] studied the laminar convection flow and heat transfer of fluid with 
constant wall temperature. Hossain and Takhar [5] analyzed the effect of radiation effects on free convection flow of a gas 
past a semi infinite flat plate. Ali et al. [1] examined the radiation effect on natural convection flow over a vertical surface 
in a gray gas. Later, Mansour [8] investigated the interaction of mixed convection with thermal radiation in laminar 
boundary layer flow over a horizontal, continuous moving sheet with suction/injection. Sattar and Kalim [10] analyzed the 
combined unsteady free convection dynamic boundary layer and thermal radiation boundary layer on a semi-infinite 
vertical plate by using the Roseland diffusion approximation. Chen [4] investigated natural convection flow over a 
permeable surface with variable wall temperature and concentration. The viscous dissipation plays an important role in the 
cooling of nuclear reactors. K. V. Chandrasekhar [3] analyzed the heat transfer of MHD flow over exponentially stretching 
sheet. Viscous dissipation effects on non-linear MHD flow in a porous medium over a stretching porous surface have been 
analyzed by S. P. Anjali Devi and B. Ganga [2]. Jha and Ajibade [6] studied the effect of viscous dissipation on natural 
convection flow between vertical parallel plates with time-periodic boundary conditions. Hassan et al [7] considered the 
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combined effect of viscous dissipation and radiation on free convection flow an inclined stretching sheet. In the present 
analysis chemical reaction was included, in addition to the problem of Hassan et al [7]. 

3. MATHEMATICAL MODEL 

A steady-state two-dimensional heat and mass transfer flow of an electrically conducting viscous incompressible 
chemically reacting fluid is considered along an isothermal permeable stretching sheet making an angle of inclination (X 
to the vertical embedded in a porous medium with heat generation/absorption. A strong magnetic field is applied along the 
y-axis direction. Here the effect of the induced magnetic field is assumed to be small enough in comparison to the applied 
magnetic field. The electrical current flowing in the fluid will develop induced magnetic field if the fluid were an electrical 
insulator, but here we have taken the fluid to be electrically conducting. Hence, only the applied magnetic field of strength 

Bq plays a role which gives rise to magnetic forces in x-direction, _ aB 0 2 in x - direction, where <J is the electrical 

u 

conductivity assumed to be directly proportional to the x- translational velocity. Two equal and opposite forces are 
introduced along the x- axis so that the sheet is stretched keeping the origin fixed as shown in Figure. 



y 



V 


Geometry of the Problem 

The fluid is assumed to be gray, absorbing emitting radiation, but non scattering medium and the Roseland 
approximation is used to describe the radiation heat flux in the energy equation. The radiative heat flux in the x-direction is 

T C 

negligible in the y-direction. The plate temperature and concentration are initially raised to w and w respectively, which 

T 

are thereafter maintained constant. The ambient temperature of the flow is 00 and the concentration of the uniform flow 

c 

is 00 . The flow and heat transfer in the presence of radiation are described as follows. 


du 3v 

- 1 - 

dx dy 


(l) 


2 2 

du dv d u 0/rr _ . O' Bn u 

u- + v—=v— + g 0 j3(T-T oo )cosa - y — 


dx dy dy : 


vu 

k 


bu z 


( 2 ) 


Impact Factor (JCC): 6.8765 


NAAS Rating: 3.11 



K. V. Chandra Sekhar & V. Manjula 


„K +V ?L = J^?L + __L 

dx dy pc p dy 2 pc p c p {dy) pc p dy 


u^- + v^ = D m ^-kr 2 (C-C co ) (4) 

dx dy dy 2 

The corresponding boundary conditions are 

u = Dx,v = v w (x),T=T W ,C = C W at y= 0, 

u = 0,T = T 00 ,C = C 00 Asy —> °°. (5) 

By using Rosseland approximation, takes the form 
4cr 3r 4 

,r=_ 3„ Sv (6) 

It is assumed that the temperature difference within the flow is sufficiently so that may be expressed as a linear 


function of temperature. 


T 4 =4TjT-3Tj 


Using (7) in (3) gives 


dx dy pc p dy 2 pc p c p \ dy J 3pc p K x dy 2 


4 . SIMILARITY ANALYSIS 

To obtain the similarity solution, the following similarity transformations are used. 


V=yA—z—,¥=^ 

v vRe v 


C' lyj ■* oo ^ VV 


Using the above transformations equations (2), (8) and (4) takes the form 

r+^f{v)f{v)+( M +^)f\v)+[i+^y^f\v) 2 - 

D y Dg) y Dg J D 


■y Re^ DO cos a 


O v (jl) + PnRe x QO + ^^f(j)) 2 + 0'(j))f(rj )=0 

D 2 D 


(p'\v)+^^f{v)(p{v)-kr 2 (p=d 


Where, 
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r= 


Gr r 


Re 


SO PVT-T_)_ = oBq 


Pr = 


X 

C p 

K 


Q = 


D 2 x 

_Qi L 

pc D 


pD 


t, k n x 2 D 

Da = — , Re x =- 


F s =~, 


Ec = 


D 2 x 2 


N = 


KK X 


(T~T~) 4(7! T; 


n 3)V Pr v 

Pn =-, Sc = — 

4 + 3N D. 


(13) 


The boundary conditions (5) transformed as 

f(rj) = Fw,f'(rj) = l,e(tj) = lp(rj) = latTj = 0 (14) 


f'{Tj) = O,0(Tj) = O,q>(rj) = OA S Tj 


—> oo 


—v 

Where Fw = ~ W 


\[dv 


Numerical Computation 

Numerical solution for the nonlinear ode (10), (11) and (12) along with boundary conditions (14) are derived 
using R - K method of fourth order through MATLAB CODE. 

5. RESULTS AND DISCUSSIONS 

For the discussion of results of flow field, the numerical computations were presented in the form of profiles of 
non dimensional velocity temperature and concentration Numerical computations have been carried out for different values 
of the suction parameter ( F w ), Magnetic field parameter (M), permeable parameter (* ) angle of inclination ( a ), heat 
source parameter (Q), radiation parameter (N), Reynolds number (Re), Eckert number (Ec), chemical reaction parameter 
( kc ), and Schmidt number ( Sc ). The results are in good agreement with That of Hasan et al [7]. The results of the 
numerical solution for the velocity, temperature and concentration profiles are depicted in Figures 1-14. 

Velocity profiles for different values of CC , k, Re, M and Fw were shown in Figure -lto 5. From Figure - 1, It was 
observed that increase of angle of inclination results in decrease in velocity, but as we move away from the sheet reverse 
trend was observed. This is due to decrease in buoyancy force. Figure - 2 depicts that increase of permeability decreases 
the velocity. But far away from the sheet the effect is not significant. Figure - 3 depicts velocity decreases with increase of 
Reynolds number. Figure - 4 shows the influence of Magnetic parameter on velocity. It can be seen that flow is slow near 
the sheet with increase of M. But later a reverse trend was observed. Figure - 5 describes effect of suction parameter on the 
flow. It is noticed that velocity decreases with increase of suction. 



Figure 1: Velocity Profiles for Different Values of Angle of Inclination 
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Figure 2: Velocity Profiles for Different Values of Porosity 



V 

Figure 3: Velocity Profiles for Different Values of Reynolds Number 



V 

Figure 4: Velocity Profiles for Different Values of Magnetic Parameter 



V 

Figure 5: Velocity Profiles for Different Values of Suction Parameter 
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Figures 6-10 shows the impact of Radiation, Heat source, Reynolds number, suction parameter and Ekert 
number It was observed that increase of Heat source, Reynolds number, Ekert number leads to increase in temperature, but 
decreases with increase of Radiation parameter, and suction. 



v 

Figure 6: Temperature Profiles for Different Values of Radiation Parameter 



v 


Figure 7: Temperature Profiles for Different Values of Heat Source 



»/ 


Figure 8: Temperature Profiles for Different Values of Reynolds Number 
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v 

Figure 9: Temperature Profiles for Different Values of Suction Parameter 



0 12 3 4 5 6 7 


Figure 10: Temperature Profiles for Different Ekert Number 

Figures 11-14 shows decrease of concentration profiles under the influence of chemical reaction parameter, 
Schmidt number, Reynolds number and suction parameter. 



Q 1 2 3 4 5 6 7 

V 

Figure 11: Concentration Profiles for Different Values of Chemical Reaction Parameter 
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Figure 12: Concentration Profiles for Different Values of Schmidt Number 
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Figure 13: Concentration Profiles for Different Values of Reynolds Number 
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Figure 14: Concentration Profiles for Different Values of Suction Parameter 

CONCLUSIONS 

The present study is to analyze the influence of different physical parameters on the MHD free convective flow 
past inclined sheet. Similarity transformations used by Hasan et al[7] for velocity components along x - direction and y- 

direction are u = Dxf \?/), v = —V Dv fifl) • But in the present discussion, the transformations given by equation (9) 

contains Reynolds number. That is why, we have shown contribution of Reynolds number on velocity, temperature and 
concentration profiles. However, the results obtained by us are in agreement with those of Hasn et al[7]. The main purpose 
of attempting this problem is due to its vast application in industries, to determine the quality of the final product. , 
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Finally, we can make the following conclusions of the present discussion. 

• Velocity decreases with increase of Reynolds number and suction parameter. But as the angle of inclination 
porosity and magnetic parameter increases velocity increases up to TJ — 2 and decreases thereafter. 

• Temperature decreases with increase of Radiation, suction parameter, but increases with increase of Heat source, 
Reynolds number and Ekert number. 

• Concentration decreases with increase of Chemical reaction parameter, Schmidt number, Reynolds number and 
suction parameter. 

Finally here, the comparative results of skin - friction ( f"( 0) ), Nusselt number ( —0 (0) ) and Sherwood 
number ( ~(p\ 0) ) are given under the influence of various parameters. 


Present Results 

Pr 

Q 

Ec 

N 
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-d'(O) 
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Results of Hasan et al 
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